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NOMI NiT AlllRl 


coolant mutiuj layon n'Ohs socUonal 
aiva 

c ii\ iimioronco ol ifu' cros^ sectional 
area 

UiMn coett 1C lent 
constant pressure spent ic heal 
major avis ot ellipse 
injection holr ili ante ter 
minor axis of ell ipse 
ompit ual constants in eq» \b) 
normal vii stance trom the wall to 
mixiiu) layer centerline 
neuiiit ot coolant mixing layer (along 
the centerline) 

Mach tutmher 

entrainment mass Liow rate 
statu pressure 
Ideal gas constant 
s i a rt i ng 1 oc a 1 1 on ol hou lul a ry 1 ay ei‘ 
t 1 ow 

local total teniperatui'e 


injeclant: toLal lemperature 
crossflow lota] temperature 
local adiabatic wall temperature 
invfscid flow velocity within the 
coolant mixing layer 
velocity within the hounuary layer 
coolant velocity 
crossflow velocity 
distance along the wall surtace 
(X « 0 at S) 

distance normal to the wall surface 
cooixiinates along the coolant iinxling 
layer centerline 

inclination angle ol the injetlion 
ho 1 e 

boundary layer thickness 
film coo 1 i ng e 1 1 ec 1 1 veiu‘s s , 




inviscid flow diivction within the 
coolant mixing layer 
dimensionless inject uiii mass flew 
rate, 

Vi scosity 

density of the coolant muimj lajcr 
coolant density 
cross/low density 
wall shear stress 






INTRODUCTION 

Interest in surface cooling, boundary layer 
control and mixing process has promoted research in 
film cooling and blowing over a flat plate. For the 
film cooling application, normal injection, tangen- 
tial slot injection, and injection through inclinea 
holes have been the subjects of theoretical and 
experimental investigations to predict the film 
cooling effectiveness, n = (To„- Taw)/(To-.- Tpr), 
and the surface heat transfer rate (refs. 1 to 5). 

In most studies of film cooling with slot or 
discrete hole injection, the local adiabatic wall 
temperature was computed from the measured heat 
transfer rate. For this process, the boundary 
layer analysis was used. However, within the re- 
gion near '•he injection, strong mixing invalidates 
the use of the heat transfer coefficient based on a 
well-developed boundary layer. 

A systematic study of coolant injection 
through a hole into a crossflow and the concept of 
an equivalent slot width to correlate all data of 
slot and hole injection was described in refer- 
ence 3. The experimental investigation of the 
effects of constant mainstream acceleration on 
effectiveness and heat transfer has shown that the 
mainstream acceleration increases the effectiveness 
downstream of a single hole (ref, 4). The interac- 
tion due to a jet entering normally into a free- 
stream has been studied with particular attention 
to the visualization of the large scale flow inter- 
action (ref. 5), and the film cooling downstream 
from a single normal jet yields a lower centerline 
effectiveness compared to an inclined jet because 
of greater mixing with freestream. In reference 6 
infrared thermography was used to measure the local 
adiabatic wall temperatures in the film cooled area 
of a plate. The test conditions involved low speed 
crossflow velocity (U« < 100 m/sec), ambient free- 
stream temperatures, an^ coolant temperatures ap- 
proximately 20 degrees below ambient. Only single 
hole injection was investigated and the coolant 
hole was inclined 30 degrees in the direction of 
the crossflow (oq = 30^*). Knowing the distribu- 
tion of wall temperature and the temperatures of 
the coolant air and the freestream, the effective- 
ness distribution, or footprint, was readily as- 
certained. 

Recently a three-dimensional finite difference 
procedure (ref. 7) has been applied in reference 8 
to the problem of predicting the flow and thermal 
fields arising from the injection of fluid in dis- 
crete jets through a wall into a crossflow. An 
example was presented to simulate a gas turbine 
blade cooled with a single row of holes aligned at 
30 degrees to the blade surface. Such an analysis 
is very complicated, time consuming on the conh 
puter, requires complex boundary conditions, and 
does not lend itself to easy parametric variations 
that reflect changes in the mixing process between 
the coolant and freestream. One of the principal 
objectives of the current research is to evaluate 
the influence of mixing on the effectiveness. 

Tfiere appears to be a need for an analytical model 
that is amenable to parametric, or modeling, 
changes. 

Reported herein is a two-dimensional integral 
method which will predict the local centerline 
adiabatic wall temperature and the coolant coverage 
area within a short dimensionless distance down- 
stream of a 30 degree inclined coolant injection 
hole. Assumptions were made about the coolant 
layer cross-section shape and the entrainment mech- 


anism Involved in the mixing of the freestream with 
tne coolant. Integral conservation equations gov- 
erning the coolant mixing layer over the surface 
were formulated and solved numerically. The re- 
sults of this analysis are compared to the low 
speed film cooling experimental results in refer- 
ence 6. 


ANALYSIS 

The problem considered in the present study is 
depicted in Fig. 1. The geometry is a flat plate, 
with a hole of diameter Dq whose axis makes an 
angle oq with the surface. Air coolant was 
Injected through this hole into an external cross- 
flow, U^. The coolant mixing layer changes its 
direction rapidly and turns towards the plate. The 
flow separates immediately after the injection 
hole. A boundary layer flow develops along the 
plate surface after the coolant layer reattaches to 
the surface. 

In this analysis, the entire mixing layer is 
first solved with inviscid flow theory. A boundary 
layer effect is then included in the reattached 
mixing layer. 

I. Inviscid Analysis 

The flow field is dominated by the mixing pro- 
cess and the pressure gradient caused by the rapid 
changing in direction in the initial stage. The 
boundary layer flow along the wall surface is neg- 
lected and the viscous effect (t^^) to the momen- 
tum conservation equation is not included. 


A. Assumptions 

(1) The flow is steady and the wall surface is 
adiabatic, 

(2) The velocity is uniform within the cross 
section normal to the axis of the mixing layer. 

(3) The coolant flow cross sectional area 
changes from a circle at the coolant exit to an 
ellipse at location bb (Fig. 1) and the ratio of 
minor axis to major axis varies linearly from 1 to 
0.25 as a function of Z/Dq. Similar elliptical 
cross sections with d/D = 0.25 exist after loca- 
tion bb. 

(4) With an adiabatic wall and low speed flow 
(M < 0.3), the local adiabatic wall temperature 
equals the local total temperature of the mixing 
layer, Taw = T. 

B, Governing Equations 


The theoretical analysis of a jet injected 
into a crossflow (ref. 9) is utilized. The energy 
conservation equation is added to solve trie local 
total temperature. 


(1) Streamwise Momentum Equation 




mU„ cos e 


With the aid of Euler's equation of the 
external flow, 

dp ii2 ■ do 

P„U„ sin 0 cos 



The streamwise momentun) equation becomes 
sin 2e ^ + AU^ ^ 


+ 2pUA pU^ s mU„, cos g (1) 

(2) Momentum Equation Normal to the C^/nterline 


first term of the right side of equation (6) and 
the extrainment due to the freestreain moving along 
the edge of the coolant layer is accounted for by 
the second term. In surface film cooling, the 
coolant layer attaches to the wall surface and a 
factor of 2 was inserted in equation (6) to con- 
sider the effective unit surface area for entrain- 
ment evaluation in the present inviscid analysis. 


The surface forces and centrifugal body forces 
are equated to the momentum flux due to entrain- 
ment. The emerging jet is treated as if it were a 
solid body in computing its contribution to the 
drag. A drag coefficient Cn is utilized for 
this purpose and the normal force is expressed as 
[Cd (pjill'i) Siri2e] D dZ 

P 

Cq “" 2 — sin 9 D ApU^ sin © (2) 

(3) Energy Conservation Equation 

With the adiabatic wall assumption, conserva- 
tion of the total enthalpy in the coolant mixing 
layer results in the following energy balance equa- 
tion: 



For constant Cp, the above equation be- 


C. Cross Sectional Area of the Coolant Mixing Layer 

Employing assumption (3), geometric relations 
for the mixing layer cross sectional area are 
derived. 

From aa to bb in Fig. 1, 


d/0 = 1 - 0.433 (Z/Oq) 


C = itD 


1 + (1 - 0.433 Z/D^y 


ll/2 


After location bb, 


(7) 


d/D « 0.25, C = 2.240, A = ir0^/16 
II. Boundary Layer Effect 


( 8 ) 


The effect of boundary layer development along 
the adiabatic wall surface after the reattadiment 
of the coolant mixing layer is considered. A two- 
dimensional boundary layer is assumed. 


PUT + pUA g + TAP 5^ + TAU ^ = niT„„ (3) 

(4) Mass Conservation Equation 

The continuity equation, with entrainment 
effect, can be written as 

(5) Equation of State 

Assuming the equation of state of an ideal gas 
is Valid in the coolant mixing layer, with the aid 
of Euler's equation, the equation of state can be 
written as 


A. Assumptions 

(1) The pressure distribution is imposed by 
the main cross flow and the streamwise pressure 
gradient is zero due to constant U„. 

(2) The wall beneath the boundary layer is 
adiabatic. 

(3) Both laminar and turbulent boundary layers 
are considered. Expressions for flat plate bound- 
ary layer parameters are valid. 

(4) The local total temperature is constant 
across the boundary layer and the local static tem- 
perature equals the local total temperature within 
the low speed mixing layer, 

B. Governing Equations 



sin 2 g 


de 


+ 


RpU dU 



Rp 


dT 

dZ 


= 0 


(5) 


where R is the ideal gas constant. 

(6) Entrainment Mass Flow Rate, m 

To solve U, p, T, A, and g from equations 
(1) to (5), the entrainment mass flow rate must be 
estimated. The following correlation similar to 
the equation used in reference 10 is employed: 

ni = p sin G 

'^00 2 “ 


+ ^|C)/2/(l £38^ sin G/U) (6) 

where E^, E2, and E3 are empirical con- 
stants. The effect of the approaching freestream 
particles to jet entrainment is modeled in the 


Principles for the development of the govern- 
ing equations are similar to section I-B. This 
technique has also been used in reference 5 for 
tangential slot injection. 

(1) Continuity Equation 

i[pU(H-«)] PU, dY^=m (9) 

(2) Z-momentum Equation 

^ [pU^(H - «)] + PU^ dY^ = + trtJ„ (10) 

(3) Energy Equation 

w ‘' y ) = ( 11 ) 


3 


(^) tquoLion of State 


tlT 


+ T 


(Jp 

Til 


a low speed coolant nilxlng layer* the local effec- 
tiveness* n « (To«. - Taw)/(TQ^. - Toe) * (To«^- T)/ 
(Tqcx>- Tgc) could be computed from these riunieri- 
(1-^) cal results. 


(i)) Entraimiient Mass Flow Rate 


RESULTS AND DISCUSSION 


For two-dimensional flow and o « 0, equa- 
tion (b) becomes 


= P.J2IU - U.„| (13) 


(b) Boumliiry Lay<2>' Parameters (ret. 11) 


For laminar boundary layer 
\ir<i 


% 

r 


a H 3 


0.332 pU 


(a) 




For turbulent boundary layer 


(14) 



(lb) 


llh Method of Solutions to the Governinci Equations 


Equations (1) to (b) were solved for the first 
derivatives of the flow variaoles as functions of 
the variables, i.e.* 


F^(p, U, I, A, g) (lo) 

where f^ = p, U, T, A, or o. 

Incliiuing the boundary layer flow effect, 
equations similar to equation (lb) can be derived 
by substituting the assumption, equations (14) and 
(lb), into equations (9) to (IJ), and resulted as 


I . In viscid Analysis 


Different values of the empirical constants, 
Cq, Ej, E 2 * and Ej have been used In 
equations (16) to find the best agreement in local 
adiabatic wall temperatures from the numerical 
solution and the experiment (ref. 6), The process 
failed to yield a unique value of Ej which will 
give satisfactory agreement for different injection 
rates. Film cooling effectiveness, based on the 
computed adiabatic wall temperature with Cq « 2, 

=* O-OU, Ej « 3U, and Ei « 0.2 and 0.3 
respectively for x « 0,3/ anu 0.b4 are compareu 
with experimental data in Figs. 2 and 3. The pre- 
sent inviscid analysis preoicts the centerline dis- 
tribution of the film cooling effectiveness over 
most of the centerline. However, there is large 
discrepancy for Z/Dq < 1. The effective cover- 
age area, computed from equations {/) and (8) are 
also shown with experimentdl data In Figs. 4 and 
b. The experimental data are based on the isothenii 
footprints for p =* 0.1. A smaller effective area 
was computed than was observed experimental 1y. 

The entrainment mass flow rates for x « 0.3/ 
and 0.b4 with prescribed empirical constants are 
compared in Fig. 6, Tfie effects of constant E^ 
and the flow direction are significant for 
Zq/Dq < 3. After Zq/Do « 3, a smaller 
entralTTinent rate was computed for x « 0.54 than 
X a 0.37 due to the velocity difference within I lie 
coolant layer. 

To investigate the present geometry assump- 
tions, the inviscid flow direction, major axis of 
the elliptical cross section, and the differences 
between the normal distance from the center of the 
analytical mixing layer to the wall and half of Liio 
minor axis are plotted for x = 0.3/ in Fig. 7, 

The flow angle, G, within the coolant layer de- 
creases rapidly from location aa to location 
bb. After location bb, a small flow angle, g, was 
computed. Very small difference, between e and 
d/2 was found. These coniputat ional results show 
tiiat the coolant layer does not separate from the 
wall surface with ag = 30 degrees and x 0.3/. 
This is consistent with the flow visualization re- 
sults in reference 12. 


= 0 .j ( p » U , H ) (17) 

where g-j p, U, or H. 

These two sets of equations were solved numer- 
ical ly as initial value problems, using the fourtn 
order Runge-Kutta method, to determine the local 
values of' p, U, T, g, A (or H). Coolant existing 
conditions were the initial conditions for equation 
(16). The inviscid flow solution from equation 
(16) at the location where g^O were used as the 
initial conditions in solving equation (17). The 
local total temperature, with the boundai*y layer 
effect, was obtained from equation (12) after the 
local density was computed from equation (17). 

With the assumption that tne local total tempera- 
ture equals the local adiabatic wall temperature in 


11, Boundary Layer Effect 

Results of the analysis with the boundary 
layer effect are also given in Figs. 2 to b. Tne 
initial conditions for the numerical solution were 
obtained from the inviscid solution corresponding 
to the location with g « 0. These locations 
corresponded to Z/D^^ 0.8 for x = 0.37 and 

2/Dq c 1.1 for X =5 0.54. When compared to the 
inviscid flow analysis, lower adiabatic wall tem- 
peratures were predicted, A laminar boundai'y layer 
assumption predicts higher effectiveness than the 
turbulent boundary layer case. This is consistent 
with less mixing in a laminar coolant layer. The 
local Reynolds number, pUX/u, within the coolant 
mixing layer was computed and found to be of the 
order of 10^ after Z/Dq ^ 2. Thus, a turbu- 
lent boundary layer assumption predicts the exist- 
ing experimental results better than the laminar 


assumption. For the present studies, expressions 
for the wholly developed boundary layer velocity 
profile, shear stress, and boundary layer thickness 
were used. To model the boundary layer effect 
downstream of a discrete hole injection, detail 
profiles of velocity, coolant concentration, and 
turbulence are required. Moreover, coolant injec- 
tion will significantly reduce the local skin fric- 
tion )n a region close to the injection. There- 
fore, the assumptions made in Section II-B may have 
contributed to the difference between the analysis 
and experiment in Figs. 2 and 3. Assuming the 
coolant mixing layer cross section is an ellipse 
with H/D = 0.25 and using the H values obtained 
from the ♦urbulent boundary layer assumption, the 
coolant coverage area was also computed and shown 
in Figs. 4 and 5. It differs slightly from the 
results of the inviscid analysis. 

III. Comparison with Existing Analysis 

Recently, analyses have been performed for 
single and multiple hole surface film cooling by 
numerical solution of the three-dimensional flow 
field (refs. 7 and 8). The analyses are effective 
in computing the average film cooling effectiveness 
far downstream of the injection hole but are time 
consuming in numerical computation. Equivalent 
technique based on tangential slot injection has 
also been proposed to compute the initial condi- 
tions for a boundary layer flow analysis due to 
inclined hole injection (ref, 3). The analysis 
presented herein is simple in computation and 
effective in predicting the centerline local adia- 
batic wall temperature and the effective coolant 
coverage area near the Injection hole. However, a 
general correlation for the entrainment mass flow 
rate, which would be valid for different injection 
rates, was, not found in this study. The assump- 
tions about the coolant mixing layer cross section 
also need further confirmation. Due to the assump- 
tion made in local total temperature and adiabatic 
wall temperature, the present method is also 
limited to a low speed coolant mixing layer. 

SUMMARY 

Assuming the local adiabatic wall temperature 
equals the local total temperature, an analysis has 
been developed to compute the centerline film cool- 
ing effectiveness and coolant coverage area beneath 
the coolant mixing layer downstream of a single 
film coolant injection hole inclined at 30 degrees 
to the crossflow. Integral conservation equations 
are formulated and solved numerically to determine 
the local adiabatic wall temperature and coolant 
spreading. Inviscid analysis, assuming elliptical 
cross section in the mixing layer, predicts the 
film cooling effectiveness within a small distance 
downstream of the injection hole. Lower adiabatic 
wall temperature (higher effectiveness) is found 
from a two dimensional analysis with a boundary 
layer effect. Smaller spanwise coolant spreading 
than the direct measurement from the experiment is 
computed from the present analytical methods. 

The present analysis is simple and its numer- 
ical results compared favorably with infrared mapp- 
ing of the adiabatic wall temperature from an 
existing film cooling experiment. However, this 
analysis is subject to further confirmation per- 
taining to the assumptions about the mixing layer 
geometry and entrainment rate. 
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Figure 2, - Centerline film cooling effectiveness, A = 0.37. 
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Figure 3. - Centerline film cooling effectiveness, X » 0. 54. 
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Figure 4. - The spanwise spread of the coolant, A => 0. 37. 
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Figure 5. - The spanwise spread of the coolant, X • 0.54. 




